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I. INTRODUCTION 

The Kramers (or Klein-Kramers) equation and the associated Smoluchowski equa- 

tion are widely discussed in the physics of Brownian motion and stochastic processes. 
Initially, the Kramers equation was introduced to describe the statistical evolution of a 
gas of Brownian particles in a potential well. It has been applied in different domains of 
physics such as colloidal suspensions, stellar systems, chemical reactions rate theory, nu- 
clear dynamics, just to mention a few [4]. The Kramers equation forms a particular class 
of Fokker-Planck equations with a constant diffusion coefficient and a linear friction 0. In 
the large time limit, the velocity distribution of the particles is close to the Maxwellian 
distribution and the evolution of the spatial density is described by the Smoluchowski equa- 
tion. The passage from Kramers to Smoluchowski equations has been the subject of many 
papers. In the early days, this passage was informal and qualitative . The derivation of 
the Smoluchowski equation from the Kramers equation is now a classical problem in many 
textbooks of statistical mechanics. In his monograph, van Kampen presents an expansion 
procedure similar to the Hilbert perturbation scheme j^. 

In this paper, we propose another derivation of the Smoluchowski equation based on the 
Chapman- Enskog expansion. Usually, the Chapman- Enskog perturbation scheme is used to 
derive the hydrodynamical equations (Euler, Navier- Stokes,...) from the kinetic Boltzmann 
equation In this context, the small parameter is played by the typical collision time 
r. The Chapman-Enskog expansion allows to determine the expression of the diffusion 
coefficients (viscosity, heat conductivity,...) in the fluid equations. In the context of the 
Kramers equation, the small parameter is played by the inverse of the friction coefficient ^. 
We can therefore expand the distribution function in powers of To first order in this 
expansion, we obtain the Smoluchowski equation. The diffusion coefficient is given by the 
Einstein relation D = kBT/^m where T is the temperature. 

We shall consider two main generalizations of this problem. First, we consider the situa- 
tion in which the ambiant medium is uniformly rotating so that the velocity which appears 
in the Kramers equation is the relative velocity with respect to the background flow. Inertial 
forces appear in the rotating frame and the expression of the Smoluchowski equation is mod- 
ified. Secondly, we consider a generalized class of Kramers equations introduced recently by 
Kaniadakis [sl, Frank p and Chavanis jlO, U] (see also an early version in Q). These 
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equations share the same properties as the ordinary Kramers equation but they are associ- 
ated with a more general form of free energy functional. This generalization encompasses 
the case of quantum statistics (fermions and bosons) and it can also find applications in 
the physics o: 



equation 
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complex media. In that context, they extend the nonlinear Fokker-Planck 
lJ | associated with the Tsallis entropy We derive here the generalized 
Smoluchowski equation from the generalized Kramers equation by using a systematic pro- 
cedure. The diffusion coefficient is given by a generalized Einstein relation which depends 
on a function implicitely determined by a second order differential equation. 

The paper is organized as follows. In Sec. m we recall the main properties of the (gen- 
eralized) Kramers and Smoluchowski equations. In Sec. IIIIl we apply the Chapman- Enskog 
expansion to the generalized Kramers equation. We account for a solid rotation of the back- 
ground medium and calculate the new terms that appear in the Smoluchowski equation due 
to rotation. In Sec. IIVI we derive a variant of the generalized Smoluchowski equation by 
using a moment method. Finally, in Sec. I3 we shortly discuss the application of these 
results to the case of systems with long-range interactions (self-gravitating systems, 2D vor- 
tices, bacterial populations,...). In that case, the Smoluchowski equation becomes non-local 
and possesses a rich mathematical and physical structure 10]. We determine the Lyapunov 
functional associated with the generalized Smoluchowski equation and interpret this func- 
tional as a free energy. In the limit of short-range interactions, this functional reduces to 
the Landau free energy and the non-local generalized Smoluchowski equation reduces to a 
generalized form of the Cahn-Hilliard equation. In Appendix we consider the case of an 
isotropic BGK operator and in Appendix El we consider a wider class of Kramers equations 
with time varying Lagrange multipliers accounting for integral constraints (energy, angular 
momentum and impulse conservation). 

II. GENERALIZED KRAMERS AND SMOLUCHOWSKI EQUATIONS 

It has been noted recently that standard kinetic equations (Boltzmann, Landau, Kramers, 
Smoluchowski,...) could be generalized so that they increase a larger class of functionals 
than the Boltzmann entropy (or Boltzmann free energy) Such generalized kinetic 

equations appear when the transition probabilities have an expression different from the one 
we would naively expect S or when the system is described by a stochastic process involving 
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a multiplicative noise depending on the density jlO|, ll6|. These generalized kinetic equations 
encompass the class of quantum kinetic equations (with exclusion or inclusion principle) 
already discussed in the literature. They can also be interpreted as "effective equations" 
trying to take into account "hidden constraints" in the physics of complex media 
The generalized Kramers equation can be written in the form jlOj : 



df df df d 

^ + V ■ — - V$ ■ — = 

dt dr dv d\ 



D 



fC"if)^ + Pfv 



(1) 



where / = /(r, v,t) is the distribution function, C is a convex function (i.e. C" > 0), 
(3 = 1/T is the inverse temperature and $(r) is an external potential. Since the temperature 
T is fixed, the above equation describes a canonical situation. We introduce the functional 



F[f] = ^J fv^dhd\ + j pMh + T J C{f)dhd^v, 



(2) 



where p = J fd^\ is the spatial density. This functional can be interpreted as a free energy 
F = E — TS, where 5* is a "generalized entropy" and E = K + W is the energy including 
a kinetic term and a potential term. When S[f] is the Boltzmann entropy, Eq. reduces 
to the ordinary Kramers equation. More generally, it is straightforward to check that E[f] 
plays the role of a Lyapunov functional satisfying F < 0. This is similar to a canonical 
version of the H-theorem. Finally, a stationary solution of Eq. is determined by the 
condition df /dt = 0, implying F = 0, and yielding 

,2 



c"(/, 



eqj 



+ <l>(r) - a 



(3) 



The equilibrium distribution feq{r, v) minimizes the free energy E[f] at fixed mass and tem- 
perature. Similarly, the generalized Kramers equation maximizes the rate of dissipation 
of free energy F at fixed mass and temperature. We refer to Ref. for a more detailed 
discussion of these results. 

For large times t ^ where ^ = DjS is the friction coefficient, the distribution function 
is close to the isotropic distribution function /(r, v,t) determined by the relation 
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C'if) = -f3 



+ A(r,t) 



(4) 



This distribution function cancels out the diffusion current in the generalized Kramers equa- 
tion (ID). The pressure P = | / fv'^d^w = p{\) and the density p = f fd^\ = p(A) are related 
to each other by a barotropic equation of state p = p{p) which is entirely specified by the 



'unction C{f). Taking the moments of the generahzed Kramers equation (^J, it is shown in 
1^ that the time evolution of the density p(r, t) is governed by the generahzed Smoluchowski 
equation 



dp 
di 



V- 



This equation was initially obtained in 



(Vp + pV$) 



(5) 



12i | from a special form of Kramers equation asso- 
ciated with the Fermi-Dirac entropy. The generalized Smoluchowski equation decreases the 
functional 

Pip') 



F[P] 



P 



P- 



(6) 



which is the simplified form of free energy ((2)) obtained by using the isotropic distribution 
function (jH) to express F[f] as a functional of p Thus, F < 0, as can be checked 

by a direct calculation. At equilibrium, F = and we obtain the condition of hydrostatic 
equilibrium 

Vp + pV$ = 0. (7) 

This relation also directly results from Eq. Q. The object of this paper is to present a 
formal derivation of the generalized Smoluchowski equation (jHI) from the generalized Kramers 
equation ((T)), using a Chapman- Enskog expansion. 



III. CHAPMAN-ENSKOG EXPANSION OF THE GENERALIZED KRAMERS 
EQUATION 



We write the generalized Kramers equation in the form 



df df df d 

— + v — -V$- — = <D 

dt dr dv dv I 



Kf)^ + Pgif)w 



(8) 



where g and h are arbitrary (positive) functions [ll| . We assume that our Brownian particles 
evolve in a background that is rotating and translating uniformly, so that w = v — u denotes 
the relative velocity of the particles with respect to the background flow with u = Oxr + U. 
The "generalized entropy" associated with this equation can be written 



S[f] 



J C{f)dVY, 



where the function C is determined by the relation 

Hf) 



C'if) 



gif) 



(9) 



(10) 



We consider the high friction hmit ^ = Dj3 +00, where (3 = 0(1). We rewrite Eq. (jS} 
in the form 

| + i/ = iQ(/). (U) 

where e = D^^ is a small parameter, L = v ■ d/dr — V$ ■ d/dv is the advective operator 
and Q{f) is the collision operator 

«<^' - 1 

Our aim is now to derive an equation for p(r, t) to second order with respect to e. For 
that purpose, we follow a Chapman- Enskog procedure and employ a formalism in the spirit 



(12) 



of 



17| . More precisely, we look for a solution / to (fTTj) of the form f = fo + e/i where /o is 



the equilibrium distribution function with the same density as /, that is, Q{fo) = and 



j /o(r,v,t)dV = p(r,t). 



Substituting this expression in Eq. (fTT|) . we get 



d_ 
di 



+ L (/o + e/0 



-Q(/o + e/i), 
e 

DQ{fo)h + Oie), 



where DQ{fo) is the linearized collision operator given by 



DQ{fo)fi 



d_ 



h'ifo)fi^ + h{fo)^ + Pg'ifo)fi 



w 



It then follows from Eq. (fT3j) that fi satisfies 



d 



DQ{fo)fi={g-^+L]fo + 0{e). 



Next, integrating Eq. (jH)) over the velocities, we obtain 

dp 



dt 



+ V 



fo^rd'\ 



(13) 

(14) 
(15) 

(16) 

(17) 

(18) 



Our aim now is to express fo and /i in terms of p, so as to deduce from (fT5|) a closed 
equation for p. First, we note that Q{fo) = is equivalent to the condition 



C'ifo) = -P 



w 



+ A(r,t) 



(19) 



This relation, supplemented by Eq. (fT^I) . completely specifies /q. 

We next determine f\. For that purpose, in view of ()17p. we need to investigate the 
invertibility properties of DQi^fo). This is the subject of the following result. 



Proposition III.l The operator DQlfo) is a non-positive self-adjoint and closed operator 
on the Hilbert space 7i = L^(]R'^, C"(/o)(i^v). The null space Af of DQ^fo) is given by 

and the range of DQlfo) is the orthogonal space to M , denoted by A/"^. 

Proof. We observe that the hnearized colhsion operator takes the simpler form 

ifo)^iC"{fo)F)^. (21) 
Consequently, for F e 7^ and G G 7i, we have 



DQ{fo)F = ^ 



{DQ{fo)F, G)n = -J gifo)^ iC"{fo)F) ■ ^ (C"(/o)G) dV, 

from which we readily deduce that DQ{fo) is non-positive and self-adjoint in Ti. Also, 
its null space is given by Eq. (|^. The closedness of DQ^fo) finally follows by classical 
arguments. □ 

We are now in a position to determine /i. Let 11 be the orthogonal projection in Ti on 
M and denote by DQ{fQ)~^ the inverse of the restriction of DQ{fQ) to A/""*". Then, applying 
successively Id — li and DQi^fo)^^ to (fT7|). we obtain 



/i = DQU,)-' {Id - n) (^^ + l)j fo + 0(6). 



(22) 



Now, noting that 



dfo 



and 



we find that 



C'ifo)^ = -/5w, (23) 



C'ifo)^ = -PVX + (3{Vufw, (24) 



C'U)t^-^t, (25) 



C'ifo 



• VA - (Vu)^w ■ u - w ■ V$^. (26) 
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In the above expressions, Vu designates the matrix [diUj) and T is the transposition, i.e. 
Aw ■ u = w ■ A^u. Since (Vu) is skew- adjoint, i.e. (Vu)^ = — Vu, we have (Vu) w • w = 
and thus, 

B 



(M_n)(| + Lj/„- ^„j^^| 



VA - (Vu) u - V$ ■ w 



(27) 



Note that the apphcation of Jrf — 11 on the first two terms of Eq. (j26p yields zero since 
they belong to N according to Proposition Ill.l. According to Eqs. (PT|) . and (j27j), the 
function /i is determined by the partial differential equation 



a(h)-^(C"(h)h) 



-(3\^X - (Vu) u - V<l> ) ■ w + 0(e) 



Introducing the pressure 



1 



and using Eqs. (|^ and ((211), fi^d after straightforward calculations that 

1. 



P 

Substituting the foregoing relation in Eq. 



:Vp = -VA , 



C'ifo) 



d_ 
dw 



9{fo)^{C"{fo)fi) 



we obtain 



P 



(Vp + pV$e//) ■ w + 0(e). 



where 



eff 



U 



is the effective potential accounting for inertial forces. 

Since Eq. fl31|) is linear in /i, we look for solutions of the form 



C"(/o)/i 



R(w)-(Vp + pV<l'e//) + 0(e) 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



Substituting this in Eq. (jSH, we find that R must satisfy the equation 



In other words, DQ{fo)iR^/C"ifo)) 



d_ 
dw 



gifo 



dRi 
dw 



-Wi. 



-w. 



/C"(/o), which yields, together with (jH^ . 



-eV • / /ivdV = V ■ 



where the matrix X = (Xij)-^^. .^^ is given by 



A'(Vp + pV$ 



eff) 



M ^ — ^ \ DQifo) 



C'ifo 



<fw = — 



Rj 



'I 

P \Cifo) 



, DQifo) 



Ri 



C'ifo 



H 



(34) 



(35) 



(36) 



In particular, X is a non-negative matrix by Proposition IIII.ll Inserting the relation ()35|) 
in Eq. p8|) . we obtain the generalized Smoluchowski equation 



| + V.(pu) = V 



A'(Vp + pV$ 



(37) 



Coming back to the equation (jHll) satisfied by R, it can be written using Eq. (fTn|) as 



-/i(/o) A^R + g'ifo)/3iw ■ V^)R = w. 



(38) 



This equation admits solutions of the form 



Rfw) 



X w, 



where x = l3w^/2 and (f){x) is a solution of the differential equation 

-/i(/o)(2x0" + 500 + ^'(/o)(0 + 2x00 = 1- 



(39) 



(40) 



We have thus established that the first order correction to the distribution function /o is 
given by 

/i = + P^*^//) ■ (41) 



and the matrix X defined in Eq. fl36|) by 



X = xh with X = 7^ 



(42) 



3p7 C^"(/o) 

Equations (j37j) . (j4(Jj) and (j42|) formally solve the problem in the general case. We shall 
now consider simplified forms of the generalized Kramers equation, corresponding to specific 
expressions of g and h hJ. First, we impose g{f) = f and h{f) = fC"{f). In that case, 
the collision operator takes the form 



fC'if)^ + /3/w 



and the differential equation pUj) becomes 

-/oC"(/o)(2x0" + 500 + + 2x0' = 1. 
Clearly, = 1 is a solution. Substituting this result in Eq. ()42j) . we obtain 



(43) 



(44) 



X 



wd^w 



3/3p J dw 
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(45) 



Therefore, for this special class of generalized Kramers equations, the Smoluchowski equation 
is 

r)n ri 1 

(46) 



| + V.(pu) = V 



-(Vp + pV$e//) 



This equation was previously derived in |lo| by working out the moments of the Kramers 
equation. The present method, in addition of being more rigorous, also provides the first 
order correction to the distribution function. We have 



/ = /o - ^T^TT-lVp + pV<l>e//) ■ w + 0{e^ 



'C"{h)p 

If we now impose g{f) = \/C"{f) and h{f) = 1, the collision operator takes the form 



d 

Qif) = IT, 



df 



p 



dv [<9v C"{f) 
and the differential equation (jiUI) becomes 



w 



-{2x(j)" + 5(f)') 



+ 2x6') = 1. 



(47) 



(48) 



(49) 



C'ifo)' 

There does not seem to be any simple solution to this equation. This implies that x is 
probably not equal to 1/^ in that case. 



IV. A MOMENT METHOD FOR THE GENERALIZED KRAMERS EQUATION 



In this section, we use a moment method to derive a variant of the Smoluchowski equation 
from the generalized Kramers equation. This approach was previously developed in 3| for 
the fermionic Kramers equation. We consider here a more general situation. We show 
in particular that the diffusion coefficient obtained by the moment method has a different 
expression from that obtained by the above Chapman-Enskog approach. To simplify the 
presentation, we consider the case u = 0, that is 



df df df d , ^ 

— + v — -V$- — = 

dt dr dv dv 



h{f)^ + P9{f)^ 



(50) 



and seek equations satisfied by the macroscopic quantities associated with / (mass, impulse 
and internal energy) defined as 



p u 



/(v) 



(51) 
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where w = v — u(r,t). Integrating Eq. (jHO)) against 1, v, and we classically get 
dtp+V- (pu) = 0, 



dtipu) + V • (puu) + V ■ P + pV$ = / Q(/)vd3v, (52) 
dtE + V ■[EvL + PvL + c^+pu-V^ = \j Qijyd^v, 

where 

E = p{— + e), P = J f{Y)wwd\, and q = y f{v)—wd\, (53) 

are the energy, the pressure tensor and the current of heat. By combining Eqs. (|52l -b) and 
fl52l -c). we can rewrite the equation for the energy as 

dtipe) + V ■ (peu) + V • q + PVu =^J Q{f)wH^w. (54) 

System ()52j) is not closed in the sense that it involves unknown quantities P{f) and q(/) 
that still depend on the distribution function /. It is well known that a maximum entropy 
principle can be used to specify / in terms of p, u and e leading to a closed set of equations 
for these quantities. Under some suitable hypotheses on the entropy functional C{f), one 
can show that the distribution function that maximizes the entropy density — J C{f )d^v 
with the prescribed conditions (jKT|) is given by 

CU) = -fci)(Y + ^(r'^))- (55) 

Then, we have 

PU*) =Ph, P=\j f*Mw'd\, and q(/,) = 0. (56) 
We note in particular that p = |pe. On the other hand, a straighforward computation gives 

Q{^)^rd'^^ = --p u, (57) 
X 



and 



with 



J QiQ'^d'^ = kp-P)J w'g{f,)d'v, (58) 



X = A c . . W3 ■ (59) 
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We thus obtain the system of hydro dynamic equations 
dtp + V ■ (pu) = 0, 

p{dtVL + u ■ Vu) = -S/p - pV$ - 4pu, (60) 
p{dte + u ■ Ve) + ■ u = i(^ - /5) / w^gU,)(P^^. 

This system of equations can be simphfied in some particular situations. In the case 
g{f) = f and h{f) = fC"{f), we find that 

X = |, j nj^9{h)d\ = 2pe. (61) 

In the case h{f) = 1 and g{f) = 1/C"{f), we obtain 

P 1 



X 



j w^g(f.)d\ = |, (62) 



where we have used a method similar to that of Eq. ()45|) to obtain the second equality. 
For the fermionic Kramers equation with h{f) = 1 and g{f) = /{rjo — /), we recover the 



equations obtained in 



g{f) = f, we have p = pT and e = |T and the hydrodynamic equations 



12l |. Finally, for the classical Kramers equation with h{f) = 1 and 

3^ 

dtp + V- (pu) = 0, 

p(9tU + u-Vu) = -V(pf)-pV<I>-epu, (63) 



{dtT + u • VT) + TV • u = 3^(T - T). 



Returning to the general case (jHUI and considering the strong friction limit e — 0, we get 
pu = 0(e), and the two last equations of (jHUI) give 

pa = x(Vp + pV$)+0(e2), and ^ = (3 + 0{e). (64) 

Inserting this expression into the continuity equation (jUUl -a). we get 

dtP = V- [xiyp + pV$)] + O(e'). (65) 

This has a similar form as the model obtained in the above section, but with new formulae for 
the quantities p and x- Now observing that the difference between the distribution functions 
/o used in the Chapman-Enskog expansion and (that maximizes the entropy) is of the 
order of e, we get p = p + 0(e). This means that, in formula (j65|) . if one replaces p by p, 
then one makes an error of the order e^ only. However, it is not clear whether X = X + 0(e) 
or not. Explicit examples are difficult to construct because of the implicit dependence of x 
in p, except for the case g{f) = f where we have x = X = V^- 

12 



V. LANGEVIN PARTICLES IN INTERACTION 



A. A non-local Smoluchowski equation 

The previous results remain valid when the force is related to the density by a relation 
of the form 

^{r,t) = j 9{r' -r)p{r',t)£r', (66) 

where 6{r' — r) is a binary potential of interaction depending only on the absolute distance 
|r' — r| between the particles. In that case, the potential energy reads 

W=^j pMh, (67) 

where the 1/2 factor guarantees that the contribution of a particle is not counted twice. With 
this form of interaction, the generalized Smoluchowski equation (j37p becomes non-local and 
takes the explicit form 



^ + u- Vp = V - |x p'{pWp + P j |^(r-r')p(r',t)rfV + pllxu |, 



(68) 



where u = ilxr + U. Equation (j68p was introduced by Chavanis JI| (see an earlier version 
in Ql). It is valid for an arbitrary equation of state p{p) and an arbitrary potential of 
interaction 9{y' — r). When the potential of interaction is long-range and attractive, this 
equation exhibits a rich variety of behaviors associated with canonical phase transitions 
and blow up phenomena. Among long-range potentials, the gravitational potential plays 
an important role. In this context, Eq. (1681) has been studied in IISI for an isothermal 

n n 

equation of state p = pT, in jl9|| for a polytropic equation of state p = Kp^ and in j20|] for a 
Fermi-Dirac equation of state. In the isothermal case, Eq. ()68|) describes a gas of Brownian 
particles in interaction (this is the canonical counterpart of a Hamiltonian system of parti- 
cles in interaction). In the non-isothermal case, it describes Langevin particles in interaction 
displaying anomalous diffusion. A formal derivation of Eq. (j68|l is proposed in [2^ start- 
ing from a generalized class of stochastic processes and using a mean-field approximation. 
Equations similar to non-local generalized Fokker-Planck equations also arise in the context 
of 2D turbulence and stellar dynamics |2^ . an d to describe the chemotaxis of bacterial 



populations 



IG 



a. AS e.p.a,ned ,n Refs. QQ, E, « ca. be used alternatively as a 



numerical algorithm to compute nonlinearly dynamically stable stationary solutions of the 



13 



Euler- Jeans equations in astrophysics. This can be particularly interesting in the case of 
rotating stars that are not spherically symmetric. 



B. The Lyapunov functional 

An additional interesting feature of the generalized Kramers and Smoluchowski equations 
when the potential $ is given by Eq. is the existence of a Lyapunov functional jlol |. 
Indeed, let us write the generalized Kramers equation in the form 



dt dv 



D 



(9v 



(69) 



where d/dt = d/ dt+L is the material derivative and J/ is the diffusion current. We introduce 
the functional 

F[f] = E[f] - TS[f] -Q L[f] - U ■ P[/], (70) 



where 



E= I f^dhd'v+^ I p^dh, (71) 



frxvdhd\, (72) 



/vd^rrf^v, (73) 

are the energy, the angular momentum and the impulse respectively. The "generalized 
entropy" S is defined by Eqs. Q and (fTUj). The functional (fTOj) can be interpreted as a 
generalized free energy. Using the skew-adjointness of Vu, it is straightforward to check 
that 

F = - J -^^d\dh < 0, (74) 

so that F[f] plays the role of a Lyapunov functional for the generalized Kramers equation. 
A stationary solution of Eq. ()69|) . defined by df /dt = 0, satisfies F = 0. Hence, it cancels 
the diffusion current = in virtue of Eq. ()74|) . It must also cancel the advection term 
Lf = 0. Using these two conditions, we can show that feq is determined by 



(75) 
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Therefore, a stationary solution of the Kramers equation extremizes the free energy ()70|) at 
fixed mass, angular velocity, linear velocity and temperature. It can be shown furthermore 
that only minima of F are linearly stable via the Kramers equation [10']. 

The Lyapunov functional associated with the generalized Smoluchowski equation (j68p 
can be obtained from Eq. (fTUI) by replacing / by its leading term /q [l0ll2J|. Thus, 



F[p] = E[fo] - TS[fo] -n Ufo] - U ■ P[/o]. 



(76) 



According to Eq. (fT^. /q can be written 



fo = F 



w 



(3[— + X(r,t) 



(77) 



where F{x) = (C) x). Using Eq. (ffTjl . the density p = f fod^^ and the pressure 
p = I J foW^cpY can be put in the form 



where 



r+oo 

G{x) = A7vV2 / F{x + t)t^/'^dt, 
Jo 



(78) 



(79) 



H(x) = ^^ F{x + t)t''^dt. 

3 Jo 

We can now express E and S as functionals of p. We follow the derivation given in 
The energy ()71|) is easily expressed in terms of hydrodynamical variables as 



(80) 



E[fo] = ljpd''^ + lj P^d'"^ +IJ P^'dh. 



(81) 



On the other hand, the generalized entropy (jH)) can be written 



d-'r 



C[F{t + p\)] t^/^dt. 



^^1^ J JO 

Integrating by parts and using C'[F{x)] = —x, we find that 



(82) 



3^53/2 J F'{t + /3A)(t + pX)e'^dt. 



(83) 



Integrating by parts one more time and using Eqs. (fTSj) . (fTUj) and (IHUj). we finally obtain 



'5[/o] = 2/5 I pdh + P I Xp dh. 



(84) 
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The angular momentum and the hnear impulse can be expressed as 

L[/o] = / prxurf^r, 



(85) 



P[/o] = J pud'r. (86) 

Collecting all these results, the free energy (f7B|) becomes 

F[p] = -J p(^\ + ^yh + ^J p<!>d'r-^ J pu'dh. (87) 

Finally, using the relation H'{x) = —G{x) obtained from Eqs. ()79|1 and ()8()|1 by a simple 
integration by parts, it is easy to check that Eq. (|7H|l implies 

(88) 

P Jo 

Hence, the free energy can be written more explicitly 

F[p] = j pf^ ^dp'dh + 1 1 p^dh - i 1 pu'dh. (89) 

This is the Lyapunov functional for the generalized Smoluchowski equation (jHEj)- Recalling 
that X = x{p) is non- negative, a straightforward calculation shows that 

- x{p)jyr < 0, (90) 
P 

where Jp is the diffusion current 

Jp = -{Vp + pV<l>eff). (91) 

At equilibrium, Jp = 0, and we obtain the condition of hydrostatic equilibrium in the 
rotating frame 

Vp + pV$e// = 0. (92) 



C. Generalized Cahn-Hilliard equation 

The free energy functional (jH^ can be written explicitly 

F[p] = j r(p)rf3r + i j j p{v)9{r' -r)p{r')dhdh' j pu'dh, (93) 
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where we have defined 

r(p)=P r^dp'. (94) 

Jo P 

Noting that 

pT"ip)=p'ip), (95) 

and taking the functional derivative of F[p], we observe that the generahzed Smoluchowski 
equation (jUHj) can be written 

| + u.Vp = V.(,pV^^). (96) 

If we now consider the case of short-range interactions, it is possible to expand the potential 

$(r,t) = y"^(r')p(r + r')rfV (97) 
in Taylor series for r' 0. Introducing the notations 

a = 4:11 6{x)x^dx and h= — / 9{x)x dx , (98) 
Jo 3 Jo 



we obtain to second order 



<l>(r,t) = ap(r,t) + ^Ap(r,t). (99) 



In that limit, the free energy takes the form 



b 



[ (Vp) 



2 



F[P] = -^ J ^^ + V{p)jd'v, (100) 

where we have set V{p) = —2T{p)/b— {a/b)p'^ + {llb)pv? . This is the usual expression of the 
Landau free energy. In general b is negative so we have to minimize the functional integral. 
On the other hand, in this case of short-range interactions, the conservative equation 
becomes 

1^ + u ■ Vp = V • l^pV (Ap - y'(p))| . (101) 

This is the Cahn-Hilliard equation which has been extensively studied in the theory of 
phase ordering kinetics. Its stationary solutions describe "domain walls". We can view 
therefore Eq. ()fi8j) as a generalization of the Cahn-Hilliard equation to the case of long-range 
interactions and arbitrary free energy functionals. This is why its physical and mathematical 
richness is so important. 
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VI. CONCLUSION 



In this paper we have derived the Smoluchowski equation from the Kramers equation by 
using a formal Chapman- Enskog expansion. We have shown furthermore that this derivation 
could be generalized to a larger class of Kramers equations where the diffusion coefficient and 
the friction term depend on the distribution function. The resulting Smoluchowski equation 
takes a nice form where the usual term Tp is replaced by a pressure p(p), as previously noted 



m 



IG 



12l |. When the potential $(r, t) is induced by the density of particles themselves, the 



generalized Smoluchowski equation becomes non-local and can generate a rich variety of 
phase transitions as discussed in Q|, and more specifically in ll^ Is^- In the limit 
of short-range interactions, it reduces to the Cahn-Hilliard equation widely studied in the 
context of phase ordering kinetics. It is interesting to note that these different equations 
are closely related to each other and that they are associated with an effective generalized 



thermodynamical formalism 
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ll| . An interest of this formalism is to unify different types 



of approaches. For example, it encompasses the case of quantum particles with exclusion or 
inclusion principle. On the other hand, our approach extends the nonlinear Fokker-Planck 
equations [l^ associated with the Tsallis entropy. Indeed, the Tsallis entropy does not play 
any special role in our formalism and most of the properties of the nonlinear Fokker-Planck 
equation can be generalized to a wider class of entropy functionals Q]. The Tsallis entropy 
forms, however, an important class of functionals associated with polytropic distributions 
and power-laws. These distributions generate a natural form of self- confinement that can 
be of interest in nonextensive systems. However, it was our interest here to show that a 
more general formalism could be developed consistently. A notion of "generalized thermo- 
dynamics", with a different presentation and a different motivation, has been developed 
independently by Kaniadakis ^, Frank and Naudts 25 1. 
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APPENDIX A: BGK OPERATOR 



In this Appendix, we derive the generahzed Smoluchowski equation from a simphfied 
form of kinetic equations. Specifically, we consider a BGK collision operator so that 

|^..§/_V*.|^^-i^, ,A1) 

ot or ov e T 

where er = er(p) is a "typical collision" time and /o is defined by 

with J fod^v = p, u = U + rixr and w = v — u. We note that contrary to the usual BGK 
operator, the velocity distribution /q is isotropic (in the rotating frame). This is consistent 
with the properties of the generalized Kramers equation (jHJ. We consider the limit e — > 0. 
Writing /i = (/ — /o) /e, we get to leading order 

fi = -r(^ + L]fo + 0{e). (A3) 



dt 

Integrating Eq. ()A1|) over the velocities, we obtain 



dt 

or, using Eq. ()A2|) . 

dt 

Now, 



^ + V ■ ( / /vc/V ) = 0, (A4) 



^ + V ■ (pu) = -eV ■ ( / fiYd^Y ] . (A5) 



(A6) 



J f,vd\=-T J v(^^ + L^/of/V + 0(e). 
The first term in the foregoing equation can be written 



9fo,3 ^ f f J3 9P (K'7\ 



To leading order, Eq. ()A5|) reduces to 

dp 



^^+V-(pu) = 0(e). (A8) 



Therefore 

dt 



v^dV = -uV ■ (pu) + O (e) . (A9) 
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The second term in Eq. ()A6|) can be written 

jm ^ I - ^ ^ / /„...*v , ,f . (MO) 

Now, 

j foViVj(fv = p5ij + pUiUj, (All) 
where p = ^ J fow'^d^\ is the pressure. After straightforward calculations, we find that 

— j hv^Vjd\ = ^ + UiV ■ (pu) + p (lixu) . . (A12) 

Combining the foregoing results, we obtain 

j f^wd\ = -riVp + pV$ + fSlxu) + 0(e). (A13) 

Substituting this in Eq. ()A5|1 and omitting terms of order e^, we obtain the generalized 
Smoluchowski equation 

dp 



g^+V.(pu) = .V 



t{Vp + pV$ 



eff) 



(A14) 



APPENDIX B: TIME DEPENDENT LAGRANGE MULTIPLIERS 



In this Appendix, we consider the generalized Kramers equation 



dt dv ^ 



Hf)^ + m9{f){^ - m^r - u(t)) 



(Bl) 



where the parameters (3, Q and U depend on time so as to satisfy at each instant t, the conser- 
vation of energy, angular momentum and impulse, and the potential $ is given by Eq. (jU^ . 
Their evolution is obtained by substituting Eq. (jBlj) in the constraints E = h = P = 0. It 
can be shown that Eq. (jBljl increases the generalized entropy at fixed mass, energy, an- 
gular momentum and impulse (H-theorem). This corresponds to a microcanonical situation. 
Equation ()B1|) can be obtained from a variational principle called Maximum Entropy Pro- 
duction Principle (MEPP). In this context, (3, O and U appear as time dependant Lagrange 
mu ltip liers associated with the integral constraints. This type of equations were proposed 
in 12] as a small-scale parametrization of the gravitational Vlasov-Poisson system in the 
context of the theory of violent relaxation. As explained in these relaxation equations 
can also be used as numerical algorithms to compute nonlinearly dynamically stable station- 
ary solutions of the Vlasov-Poisson system. This would be particularly interesting in the 
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case of rotating stellar systems that are non spherically symmetric. They can have similar 
applications in other domains of physics with different types of interaction. 

Considering the limit ^ — > +00, the Chapman-Enskog method of Sec. Illll can be im- 
plemented with only slight modifications. Due to the time dependence of the Lagrange 
multipliers, Eq. is replaced by 

C"(/o)§^ = -af +,aw.u-;3(!^ + A), (B2) 

where /o is still given by Eq. (fT^ and u = U+Oxr. Consequently, introducing the notation 



and using the identity 



resulting from Eq. (f^ . Eq. ^T7\i becomes 
,«_n)(| + L)/„.-^(vA-u-(Vu,u-V*).w-^(^-i^).(B5) 



Using Eqs. ()B5|) . we find that Eq. (|2g|) is replaced by 

/i = -P DQifo)-' (y^F^ (VA - u - (Vu) u - V<l>) ■ w 

Since /o(w) is spherically symmetric, only the terms proportional to w in the right-hand 
side of Eq. ()B6|) will give a non- zero contribution in Eq. ()18p. Therefore, Eq. ()37p remains 
valid with the substitution V$ — > V$ + u. We thus obtain the generalized Smoluchowski 
equation 

dp r 1 

^ + V-(pu) = V- x(Vp + pV$e// + pu)J. (B7) 

Note that the last term was forgotten in ^L^. The time dependent Lagrange multipliers 
Pit), Q,(t) and U(t) are determined by replacing / by its leading term /o in the conservation 
of energy, angular momentum and impulse. Recalling (jHT|) . (jH3j) and (jHBj) . we obtain 

^ j p{p)d'^r + ^ y P^d^r + ^ y pu^dh = const. , 

prxud'^r = const., 

pud^r = const. , 
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where p{p) is given by (fTHj) . (f7^ and 
temperature has been studied in [3, lldl l26 



. This type of equations with time dependant 
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